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We investigate the probability density of rescaled sum of iterates of sine-circle map within quasi- 
periodic route to chaos. When the dynamical system is strongly mixing (i.e., ergodic), standard 
Central Limit Theorem (CLT) is expected to be valid, but at the edge of chaos where iterates have 
strong correlations, the standard CLT is not necessarily to be valid anymore. We discuss here the 
main characteristics of the probability densities for the sums of iterates of deterministic dynamical 
systems which exhibit quasi-periodic route to chaos. At the golden-mean onset of chaos for the sine- 
circle map, we numerically verify that the probability density appears to converge to a g-Gaussian 
with g < 1 as the golden mean value is approached. 
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I. INTRODUCTION 

The standard Central Limit Theorem (CLT) asserts that, under appropriate conditions, the probability density of 
the sum of large number of independent identically distributed {iid) random variables will be approximately normal (or 
Gaussian) [1]. The first complete proof was given in 1901 by Liapunov who worked analitically with characteristic 
functions. Lindeberg gave a complete proof of the CLT under more general conditions than Liapunov in his 1922 
paper [3]. The CLT explains why many stochastic processes are of relevance in physics, chemistry, biology, economics, 
etc. One of the fundamental importance of CLT for statistical mechanics is that if a CLT is valid for the driving 
forces in a many-body system, it is easy to proceed to the formalism of statistical mechanics via the Langevin and 
Fokker-Planck approaches. 

It is also well-known that there are CLTs for the deterministic dynamical systems. Although deterministic dynamical 
systems have deterministic algorithms and therefore never be completely independent, if the assumption of iid is 
replaced by the weaker property that the dynamical system is sufficiently strongly mixing, then various CLTs for 
these systems are also shown to be valid 'i'-tq. Here, "sufficiently strongly mixing" implies ergodicity and just means 
asymptotic statistical independence for large time differences. In this case the standard CLT is valid for deterministic 
dynamical systems. If the iterates of deterministic dynamical systems have strong correlations like at the edge of 
chaos where the system is only weakly chaotic, then the standard CLT is not valid. In fact, many complex systems 
in nature, such as financial data earthquakes [loj . long-range interacting many body systems |Tl| . exhibit strong 
correlations among appropriate random variables. At least for a class of certain strong correlations (referred to as 
q- independence), it has been analytically shown that the attractors are g-Gaussians leading to a q-genaralization of 
the standard CLT [T^ - [T5| . g-Gaussian distributions, defined as, 

F(x) = [(l-(l-q)/3a;2]V(i-«), (1) 

(where /3 controls the width of the distribution) are the distributions that optimize the nonadditive entropy Sq (defined 
to be Sq = {\ — X^jP?) I {l~ l))j which nonextensivc statistical mechanics is based [lB,[l3- As g — > 1, g-Gaussians 
recover the Gaussian distribution. 

In recent years, many dissipative deterministic dynamical systems like cubic map, logistic map and logaritmic map 
are numerically investigated in terms of their central limit behaviour jT8j-|22| . The common property for all of these 
examples is that they exhibit period doubling route to chaos. Another dissipative but high dimensional system known 
as Kuramoto model has also been studied within the same context p3| . Finally, a very recent work has also appeared 
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as a first attempt of a similar analysis of conservative maps [2^ . In this work, we focus our attention to another one- 
dimensional dissipative dynamical system, known as sine-circle map (25l . [26| , which differs from the above-mentioned 
maps since it exhibits quasi-periodic route to chaos. We are mainly interested in two important questions: (i) what 
are the typical probability distributions in the vicinity of the quasi-periodic edge of chaos? (ii) what is the relationship 
between the degree of quasi-pcriodicity of trajectories and the complete shape of appropriate limit distribution? 
Let us start by introducing the sine-circle map 



Ot+i = 9t + n-—sm {2TT9t) mod (1) , (2) 

ZTT 

where < < 1 is a point on a circle and parameter K (with if > 0) is a measure of the strength of the nonlinearity. 
It describes dynamical systems possessing a natural frequency wi which are driven by an external force of frequency 
LU2 (fi = LOxjiOi is the bare winding number or frequence-ratio parameter) and belongs to the same universality class 
of the forced Rayleigh-Benard convection j^^l- Winding number for this map is defined to be the limit of the ratio 



11^(^1^, (3) 

t^oo t 

where (9t — 9q) is angular distance traveled after t iterations of the map function. The map is monotonic and invertible 
(nonmonotonic and noninvertable) ioi K < 1 {K > 1) and develops a cubic inflexion point at 6* = for if = 1. The 
winding number 11^(57) can be numerically computed from Eq. ([3]) forming a "Devil's Staircase" shown in Fig. 1. Once 
mode is locked, it is noted that W does not change for a substantial range of f2. If f2 belongs to a constant interval of 
the Devil's staircase (wi /w2 is commensurate) , then W is rational number (is also called mod- locked) and the behavior 
of the system is periodic, otherwise, i.e., if fl does not belong to a constant interval (cji/a;2 is incommensurate), then 
W is irrational number and the behavior of the system is quasi-periodic. 

At the onset of chaos where if = 1, a set of zero measure and universal scaling dynamics is produced at special 
irrational dressed winding numbers. They can be approximated by a sequence of truncated continued fractions. The 
most interesting and well-studied case is the sequence of rational approximants to Wqm = (V^— l)/2 which is called 
the golden mean and this has the form of an infinite continued-fraction 



W= ,^ 1 • (4) 

If fraction lines stop at n denominator, it is the nth order approximation to the golden mean. It is also easy to see 
that the nth approximation to the golden mean is given by Wn = {Fn/ F„+i), where Fn is nth Fibonacci number and 
Wn is the nth convergent for the golden mean. This states that the sequence of rational numbers Wn converges to 
the irrational number Wgm as n — > cxo, yielding the frequency-ratio parameter to approach its limiting value ftoo- 

This map has already been studied in at the edge of chaos, where the sensitivity to initial condition behavior 
has numerically analyzed for the first time and its connection with the nonextensive statistical mechanics has been 
established. Then, the numerical results given in (28| have been analytically proven by Hernandez- Saldana and 
Robledo |29|] . In the remainder of this work, we focus on the probability density of rescaled sum of iterates of the 
sine-circle map in the strongly chaotic regime and at the quasi-periodic edge of chaos and discuss the main features 
of the probability distribution as the quasi-periodic edge of chaos is approached. 



II. PROBABILITY DENSITIES OF THE SUMS OF ITERATES OF THE SINE-CIRCLE MAP 



A. Strongly chaotic case 

Generally, one can consider the sine-circle map given in Eq. Q and investigate the behavior of this dynamical system 
where it is strongly chaotic. When this is the case, implying asymptotic statistical independence, then correlations 
among distributed random variables asymptotically converge to zero, and thus the standard CLT is expected to 
be valid for this case. In order to numerically illustrate this, we can consider the sine-circle map at if = 5 and 

= 0.606661063469 where the map is strongly chaotic and check whether the probability distribution of 
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FIG. 1: Devil's staircase which shows the mode-locking structure of the sine-circle map. Some representative rational numbers 
for which the winding number is locked are given on the curve. 
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FIG. 2: (Color online) (a) Probability distribution function of the sine-circle map for A" = 5 and Q. = 0.606661063469. Black 
dashed line corresponds to the Gaussian P(y) — Z'^"" ^/V2tto^. (b) Logarithmic plot of the same data where a linear 
tendency, which corroborates the Gaussian behavior, is appreciated. 



becomes Gaussian for N oo, after appropriately scaled and centered, regarding the initial value as a random 
variable. Here, the average (...) is calculated as time average. The result is given in Fig. 2 where a clear convergence 
to a Gaussian is evident as expected. 



B. Approaching the quasi-periodic edge of chaos 

We are now at the position to investigate the behavior of the sine-circle map in the vicinity of the quasi-periodic 
edge of chaos for which the standard CLT is not expected to be valid due to strong correlations among the iterates of 
the map. At the critical point the system is not ergodic and strongly mixing anymore. To numerically calculate the 
averages (...) in Eq. ([5|), due to lack of ergodicity, it is necessary to take the average over not only a large number of 
A'^ iterations but also a large number of M randomly chosen initial values 
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TABLE I: Numerical values of used in the simulations and the associated values of n, W„ and N*. 



n 






11 


0.606664795117167 


144^ 


13 


0.606661528216938 


377^ 


15 


0.606661121349764 


987^ 


oo 


0.606661063469... 


<x 



M N 

j=i i=i 

Although this condition is of no importance at strongly chaotic regime because of ergodicity, this is certainly of great 
importance at the edge of chaos due to non-ergodicity. 

At the exact quasi-periodic edge of chaos, we know that K ~ 1 and the frequency-ratio parameter should be taken 
as an irrational number, like fioo, corresponding to Wgm- Unfortunately, in a numerical experiment, this parameter 
can never be obtained with infinite precision (which means that the winding number is exactly equal to the golden 
mean value) but can be approached with increasing precision (namely, the winding number is getting closer to the 
exact golden mean value as larger Fibonacci numbers are used). In this case, since the golden mean is approached with 
the ratio of the Fibonacci numbers Wn = (F„/F„_|_i), the trajectory has the period Fn+i- It means that trajectory 
returns to its starting point on the circle. Therefore, there exist Fn+i quasi-periodic bands of the attractor within its 
self-similar structure. As a result, the critical behavior of the map can never be numerically studied at the exact chaos 
threshold where the frequency-ratio parameter is expected to be completely irrational corresponding to the case of W 
approaching to the golden mean with infinite precision, but can only be analyzed in the vicinity of the critical point, 
systematically getting closer to it by increasing the finite precision of the golden mean adjusting the frequency-ratio 
parameter. This procedure is exactly the same as the one used in p^ . [20| for the same analysis of the logistic map. 
Similar to the case of the logistic map, one can write down a scaling relation [2^ 



\n,,{K)-n^{K)\^j^ , (7) 

where \6\ = 2.83362... is the Feigenbaum constant, ilao is the value whenever Wgm is achieved. In practice we localize 
the f2„ values from the scaling relation ([7]) and approach the critical point by increasing n values. Numerical values 
used throughout this work have been given in Table I. 

As a result of above-mentioned procedure, we are not exactly at the chaos threshold and therefore the sum of 
iterates ^iH essentially approach to a fixed value v = Fn+i{9) plus a small correction Avi which describes 

the small fluctuations of the position of the (_F'„+i)th iterate within the quasi-periodic band. Therefore, one can write 

2/1= ^{e,^{e))^Av, , (8) 

i=l 

and 

2F„ + i 

2/2 = J2 - (^)) = . (9) 

i=_F„+i + l 

If we continue this iteration Fn+i times respectively, we obtain Fn^i strongly correlated random variables 
(2/I1 2/2j ■■•5 2/-F„+i) = {Avi, Av2, Avp^^^-^) at last. Every new fluctuation Avj is correlated with the others and 
correlations among these random variables decay very slowly if we are close to the irrational frequency-ratio param- 
eter. Also, we use the total of these strongly correlated random variables to obtain appropriate limit distribution 
function in the vicinity of the quasi-periodic edge of chaos 
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FIG. 3: (Color online) (a) Probability distributions of the cases f2ii, Q.13 and fiis with appropriate A''* values, (b) q-logarithmic 
plot of the same data for the case 57i5 which is the closest to the critical point that we can numerically obtain in our simulations. 



j—1 i—1 i—1 

This clearly implies that, in order to see the convergence to the limit distribution, the appropriate number of iterations 
N* must be taken as 

N* - (F„+i)2 . (11) 

Our main results are given in Figs. 3 and 4. We plot the probability density function in Fig. 3a for f2„ values 
approaching ^l^o for the first three cases given in Table I with n = 11, n = 13 and n = 15 with appropriate N* values 
obtained from the scaling relation. It is clearly seen from Fig. 3a that the curves have a tendency to converge to a 
limit distribution as the critical point is approached. It seems that the limit distribution can be well approximated 
by a g-Gaussian with q = 0.925. When n is small (more distant case to the critical point), the data converges to 
this g-Gaussian distribution at the more central regions but deviates in the tails. As n increases (getting closer to 
the critical point), this convergence develops towards the tails. It is clear that for a perfect convergence, one needs 
to achieve n — > 00 and N 00 limits. In order to further strengthen these results, we also plot in Fig. 3b the 
q-logarithm [defined to be log^(x) = (a;^~* — !)/(! — q)] of the case which is the closest to the critical point that we 
can numerically obtain. A satisfactory agreement between the data and (?-Gaussian can easily be appreciated for the 
case with the largest n (better convergence in the tails could have been achieved if cases that are closer to the critical 
point would have been numerically attainable). In Fig. 4 we also give the same demonstration of Fig. 3 for the closest 
case to the critical point where we have used fioo — = 0.606661063470 with three different N values which are 
smaller than N* since the value of N* is beyond the present computational limits. As expected, the limit distribution 
would converge to the appropriate g-Gaussian for the entire region whenever N = N*. Unfortunately, for the present 
case, since we are very close to the critical point, the appropriate N* value is far beyond our computational limits. 
Nevertheless, as N values increase, the convergence to the appropriate g-Gaussian (which is with g = 0.90), starting 
from the central part, can easily be seen. For the largest N value that we can numerically attain, almost the entire 
region (the central part and the tails) is well approached to the g-Gaussian distribution. 

The observed behavior obtained here for the sine-circle map appears to be exactly similar to that of the logistic 
map in the sense that the numerical results agree better and better with the appropriate g-Gaussian as the critical 
point is approached although the route to chaos here is via quasi-periodicity whereas it is via period doublings in the 
logistic map. The only difference is that the g values obtained here are below unity, whereas for the logistic map they 
are always above unity. 



o 



' 1 


1 






: M=5.10'' 






- 


/ 






(a) : 


/ 

/ 




f \ 

\ \ 


: 

- 


i 




i| 










1 


— - q-Gaussian; q=0.90 ; (3=2.95 


_ \ 




2 -1 


yP(0) 


1 







1 


- - q-Gaussian; q=0.9() ; P=2.95 


N 

- 




(b) 


M=5.10' 

1 


1 


N 

1 



[yP(0)] 



FIG. 4: (Color online) (a) Probability distributions of the case Qc with increasing A'^ values with < A'^* (for this value of 
n the appropriate A'^* value is far beyond the present computational limits), (b) 5- logarithmic plot of the same data for the 
largest possible TV that we can numerically obtain in our simulations. 



III. CONCLUSIONS 



Main results obtained in this work can be summarized as follovifs: (i) The sine-circle map, which belongs to the 
same universality class of the forced Rayleigh-Benard convection, is investigated and it is shown that its appropriate 
probability distribution approaching the quasi-periodic edge of chaos appears to be q-Gaussian. Up to now, in the 
literature, the only chaos route investigated within this procedure was the period doubling route where the q values 
seem to be always larger than unity. Since the quasi-periodic route to chaos is different from the period-doubling 
route to chaos, our analysis performed here can be thought as a step forward towards the inclusion of another route 
to chaos (with q < 1) to the phenomenon observed so far for the period-doubling route to chaos (with q > 1). In this 
sense, our results broaden the applicability region of the arguments argued in [20l | and provide further evidence 
to the proposition that, for certain classes of strongly correlated random variables, the rescaled sum of the iterates of 
such dynamical systems approaches a q-Gaussian. 

(ii) Relationship between quasi-periodicity degree of trajectories and complete shape of appropriate probability 
distribution is investigated. It is seen that the probability distribution converges to a g-Gaussian as irrationality 
degree of frequency-ratio parameter increases, which also means that the winding number approaches the golden 
mean value (in other words, this is to say as the critical point is approached). 
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